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Abstract 

A quantum field-theoretical model, which describes spatially non-homogeneous repul- 
sive Bose gas in an external harmonic potential is considered. Two-point thermal 
correlation functions of the Bose gas are calculated in the framework of the functional 
integration approach. Successive integration over the "high-energy" functional vari- 
ables first and then over the "low-energy" ones is used. The effective action functional 
for the low-energy variables is obtained in one loop approximation. The functional 
integral representations for the correlation functions are estimated by means of the sta- 
tionary phase approximation. A power-law asymptotical behaviour of the correlators 
of the one-dimensional Bose gas is demonstrated in the limit, when the temperature is 
going to zero, while the volume occupied by the non-homogeneous Bose gas infinitely 
increases. The power-law behaviour is governed by the critical exponent dependent on 
the spatial arguments. 
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1 Introduction 



Experimental realization of Bose condensation in vapours of alcali metals confined in the 
magneto-optical traps stimulated a considerable interest to the theory of the Bose gas [1,2]. 
In particular, Bose condensation in the systems, which are effectively two-dimensional or 
quasi one-dimensional, became a subject of experimental and theoretical investigations. For 
more details one should be referred to [1,2]. The field models, which describe the Bose 
particles with delta-like interparticle coupling confined by an external harmonic potential, 
provide a reliable background for theoretical description of experimental situations [1,2]. 
For a translationally invariant case, the field models in question correspond to a quantum 
nonlinear Schrodinger equation which allows to obtain closed expressions for the correlation 
functions in the one-dimensional case [3]. 

Some of the results of the papers [4-6] devoted to the correlation functions of the weakly 
repulsive Bose gas confined by an external harmonic potential arc reported below. Since 
there are no exact solutions in the case of an external potential, the functional integration 
approach (see [7-13] as a hst, though incomplete, of appropriate refs.) is used in [4-6] for 
investigation of the two-point thermal correlation functions. It will be demonstrated below 
that the presence of the external potential results in a modification of the asymptotical 
behaviour of the correlation functions in comparison to a translationally invariant case. 

The paper is organized as follows. Section 1 has an introductory character. A description 
of the one-dimensional model of non-relativistic Bose field in question, as well as a summary 
of the functional integration approach, are given in Section 2. Successive integration over the 
high-energy over-condensate excitations first and then over the variables, which correspond 
to the low-energy quasi-particlcs, is used in the given paper for a derivation of one loop 
effective action for the low excited quasi-particles. The method of stationary phase is used 
in Section 3 for approximate investigation of the functional integrals, which express the two- 
point thermal correlation functions. Specifically, the asymptotical approach to estimation 
of the correlators, which is discussed in the present paper, was proposed in [14]. It is clear 
after [4-6] that the method [14] admits a generalization for the spatially non-homogeneous 
Bose gas in the external potential as well. The asymptotics of the two-point correlation 
functions of the non-homogeneous one-dimensional Bose gas are obtained in Section 4. A 
short discussion in Section 5 concludes the paper. 

2 The effective action and the Thomas— Fermi approx- 
imation 

2.1 The partition function 

Let us consider one-dimensional repulsive Bose gas on the real axis R 9 x confined by an 
external potential V{x). We represent its partition function Z in the form of the functional 
integral [7-13]: 




(1) 
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where S[ip,ip] is the action functional: 

S[^,^] ^ j\r j dx[i,{x,T)[^-ny{x,T)- 

"2 ^(^' '^)^(^' r)ip{x, r)| , 

and Vt/jV-ip is the functional integration measure. Other notations in (1), (2) are: H is the 
"single-particle" Hamiltonian, 

m is the mass of the Bose particles, /i is the chemical potential, g is the coupling constant 
corresponding to the weak repulsion (i.e., g > 0), and the external confining potential is 
V{x) = ^mQ?x'^. The domain of the functional integration in (1) is given by the space of 
the complex- valued functions r), ip{x,T) depending on a; G M and r G [0,/3]. With 
regard to x, the functions '?/'(a;,r), ^(x, r) belong to the space of quadratically integrable 
functions -L2(M), while they are finite and periodic with the period (3 — (kBT)'^ with regard 
to the imaginary time r (/cs is the Boltzmann constant, and T is an absolute temperature). 

At sufficiently low temperatures each of the variables ip{x,T), ip{x,T) is given by two 
constituents: 

ljj{x,T) ^ ljjo{x,T) +1pe{x,T) , 'ljj{x, t) ^ 'lpo{x, t) + Tpeix, t) , (4) 

where 'ipo{x,T), ipo{x,T) correspond to quasi- condensate (a true Bose condensate does not 
exist in one-dimensional systems [11]), while 'ipe{x, t), ipeix, t) correspond to the high-energy 
thermal (i.e., over-condensate) excitations. In the exactly solvable case, the existence of 
the quasi-condensate implies that a non-trivial ground state exists [3]. Let us require the 
variables (4) to be orthogonal in the following sense: 

J 7jjo{x,T)'ljje{x,T)dx ^ J o{x , t)^^ e{x , t) dx = . 

Then, the integration measure VipVip is replaced by the measure 'Dipo'^i'oT^'^eDi'e- 

To investigate the functional integral (1), we shall perform a successive integration: first, 
we shall integrate over the high-energy constituents given by (4), and then over the low- 
energy ones [8,11]. At a second step, it is preferable to pass to new functional variables, which 
describe an observable "low-energy" physics in a more adequate way. After the substitution 
of (4) into the action (2) we take into account in S only the terms up to quadratic in 'i/'e, ipf.- 
This means an approximation, in which the over-condensate quasi-particles do not couple 
with each other. In this case, it is possible to integrate out the thermal fluctuations ^e, '4^e 
in a closed form and thus to arrive to an effective action functional 5'eff['0o, ^o]- It depends 
only on the quasi-condensate variables V'o; V'o- 

SA^o. i^o] = In y e^I^o+'^-^o+^^lPVe^PV^e , (5) 

where the tilde in S implies that "self- coupling" of the fields ifjei ipe is excluded. With respect 
of (5), the partition function of the model Z {1) takes an approximate form: 



/ 



e^efilV-o.^olp^^p^^ . (6) 
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Let us consider the derivation of the effective action Ses['il>o:'4>o\ (5) in more details. The 
sphtting (4) allows to derive -S'esfV'o, V'o] in the framework of the field-theoretical approach 
of loop expansion [15]. We substitute (4) into the initial action S'['0,'0] (2) and then go over 
from S to the action S, which is given by three terms: 

+ + (7) 

In (7), iScond is the action functional of the condensate quasi-particles, which corresponds to 
a tree approximation [15]: 

'S'cond [^o,^o] = / dr / dx\'i/jo{x,T)K+'ilJo{x,T)- 

Jo J (8) 
9 " - 1 

-- ^o{x, T)7jjo{x, T)ljjo{x, T)'tjjo{x, t) |. 

At the chosen approximation the action for the over- condensate excitations 5'free takes the 
form: ^ 

5free[^e,4] = dr J da;(V^e, ^e)^-^ , (9) 

where is the matrix-differential operator, 

G-^ = Go^-E. (10) 

In (10) we defined: 

^ d 

where K± arc the differential operators, Kj^ = ±— Ti,. and Ti is the single-particle 

or 

Hamiltonian (3). Eventually, 5'int describes a coupling of the quasi-condensate to the over- 
condensate excitations: 

Sint['tpo,'^o,''Pe,'tpe] = f dr [ dx\ ^e{x, t)[K+ - g'lp„'ljjo\'lp„{x, t) + 

Jo J ^ (12) 

+1pe{x,T) [K_ - gipo'lpo\'lpo{x,T)y 

It is appropriate to apply the stationary phase method to the functional integral (6). To 
this end, let us choose 'ipo, ipo as the stationarity points of the functional Scond (8), which are 
defined by the cxtrcmum condition 5(S'cond[V'o, '0o]) — 0. The corresponding equations look 
like the Gross-Pitaevskii-type equations [1]: 



(13) 



The contribution of the action functional Si^t (12) drops out from (7), provided ■0o, ■^o are 
solutions of equations (13). Therefore the dynamics of ■^g, ■0e is described, in the leading 
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approximation, by the action ^free (9). The latter depends on -00, ■^o non-trivially through 
the matrix of the self-energy parts E, which enters into (10). 

The Thomas-Fermi approximation is essentially used in the present paper in order to 
determine the stationarity points -00, "00 ■ This approximation consists in neglect of the 

kinetic term — 7— in equations (13) [1,2]. The Thomas-Fermi approximation is valid 

2m ox^ 

for the systems containing a sufficiently large number of particles in the magneto-optical 
traps [1,2]. The following condensate r-independent solution can be obtained: 

Mo = Ptf{x-^) = -g{l^- V{x))Q{ii - V{x)) , (14) 

where © is the Heavyside function. Now the integration in (6) with respect to '0e, "06 is 
Gaussian. This leads to the one loop effective action in terms of the variables ■0o, i^o'- 

Sesii^o, ^o] = ScondbPo, ipo] " | In Dct (p'^^ . (15) 

Here is the matrix operator (10) and ipo-, "00 have a sense of the new variables governed 
by the action (15). 

In order to assign a meaning to the final expression for the effective action (15), it is 
necessary to regularize the determinant Det (^G~^^. In our case, the operator is already 

written as 2 x 2-matrix Dyson equation (10), where the entries of E['0o, i/So] (11) play the role 
of the normal (En = E22) and anomalous (E12, E21) self-energy parts. The Dyson equation 
(10) defines the matrix G, where the entries have a meaning of the Green functions of the 
fields ipe, ipe- The matrix G arises as a formal inverse of the operator G~^. 
It is appropriate to represent G~^ as follows: 

G-' = - E = - (e - 2gpTF{x; //)/) , (16) 
where / is the unit matrix of the size 2x2, and the matrix is defined as 



K_-2gpTpix;fi) \0 /C_ 



(17) 



Here ptf{x; fi) is the solution (14) and equation (16) implies that we simply added and 
subtracted 2gpTF{x', p) on the principle diagonal of the matrix operator G^^. A formal 
inverse of can be found from the following equation, which defines the Green functions 
G±: 



(18) 



/c_y \o g 

Using the relation In Det = Tr In, one gets: 



~ In Det (g-') = -1 IV In (/ - ^(S - 2^ptf(x; -\ In Det {^^ ^ ^ 



The first term in right-hand side of (18) is free from divergencies. Let us consider the 
determinant of the matrix-differential operator in right-hand side of (18). Let us denote 
the eigenvalues of the operators K,± as iia;^ — A„, where ujb are the bosonic Matsubara 
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frequencies and A„ are the energy levels labeled by the multi-index n [5,6]. Then, we 
calculate [5,6]: 



4lnDet^^- 



2/3 \0 IC 

where has a sense of the free energy of an ideal gas of the over-condensate excitations. 
Then, in the leading order in g, one gets: 

-^InDet (d-'^ ^ -f3Kc{p) + 

+g j dr j dx{g+{x,T;x,T) + g-{x,T;x,T)){'ijjoipo- Ptf{x;ij.)) = 

= -PFM - 2^ dr j dx Pnc{x)^o^o ■ (19) 

Here Fnc is the free energy of the non-ideal gas of the over-condensate quasi-particles. The 
density of the over-condensate quasi-particles is Pnc{x) = —Q±{x,t;x,t), and it depends 
only on the spatial coordinate x. At very low temperatures and sufficiently far from the 
boundary of the domain occupied by the condensate, the quantity Pnc(^) can approximately 
be replaced by Pnc(O), since Q±{x,t;x,t) is almost constant over a considerable part of the 
condensate [16]. 

It is appropriate to write the one loop effective action obtained in terms of new inde- 
pendent real- valued variables of the functional integration. Namely, in terms of the density 
p{x, t) and the phase (p{x, r) of the field ipo{x, t): 



V^„(x,r) = v/;^(^c^'^(^'^) , i,o{x,T) = ^r^{^)e-'^^^'^^ . (20) 
In terms of these variables, the effective action takes the form [5,6]: 



5'eff[p, ^] = -(3Fr,c{p) + i ^ dr J dxi^pdr^p + ^ dx{pdxip)^ + 

f^r J dx!^^[^pdl^p- pidx^r) +{A-V)p-^-p'y 



(21) 



where A — p — 2gpj^c{0) is the renormalized chemical potential. Here and below we denote 
the partial derivatives of the first order over r and x as dr and dx, respectively, whereas the 
partial derivatives of the second order — as and d^. The model in question in the present 
paper is spatially one-dimensional, and a possible multi-valuedness of the angle variable (p 
is left aside. 

We shall consider SeS (21) as the one loop effective action, where the thermal correc- 
tions over the "classical" background are taken into account. The "classical" background 
corresponds to the solution (14). It should be noticed that our derivation of the effective 
action can formally be used for two and three dimensions also [4]. Notice that equation (21) 
remains correct V — also. 

2.2 The excitation spectrum 

Let us determine the spectrum of the low-energy quasi-particles. Now we are applying 
the stationary phase approximation to the integral (6), where the effective action is given 
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by (21), while the measure is VpVip. The corresponding stationarity point is given by 
the extremum condition S(^Sef[\p,(f\) — 0, which is equivalent to the couple of the Gross- 
Pit aevskii equations: 

^drV^ + ^(^ 91Vp - idM') + A - V{x) -gp = 0, 

2m VVp J (22) 

-idrP H {pdx(p) = . 

m 

We use the Thomas-Fermi approximation and drop out the term (<9|y^)/^/p in the first 
equation in (22). Solution with drP = = appears, provided the velocity field v = 
m~^dx(p is taken equal to zero in (22). In this case, equations (22) lead to the density of the 
condensate: 

A~.. A/ x'^ \ ^ f x'^ 



Pxp(x).-?xpW = -(l-^je(l-^j. (23) 

Explicit form of the external potential V{x) = ^mfl'^x^ is taken into account in the ex- 
pression (23). The form of the solution (23) means that the quasi-condensate occupies the 
domain \x\ < Rc at zero temperature. The length R^ defines the boundary of this domain, 
2A 

i?^ = — — (in three dimensional space, this would correspond to a spherical distribution 

of the condensate). In the homogeneous case given by the limit 1/Rc 0, the Thomas- 
Fermi solution ptf{x) is transformed into the density Ptf(O) = A/g, which coincides with 
the density of the homogeneous Bose gas [8] . 

According to the initial splitting (4), we suppose that thermal fluctuations in vicinity of 
the stationarity point (23) are small, and therefore an analogous splitting can be written for 
the condensate density also: 

Po{x, t) = Ptf{x) + TToix, t) , (24) 

where Po{x,t) is a speciflc solution of (22). We linearize equations (22) in a vicinity of the 
equilibrium solution po — PTpi^), (p — const. Eliminating the phase (p and dropping out the 
terms proportional to /i^, we go over from (22) to the Stringari thermal equation [17]: 

^a^o + a.(^(l-^)a.7ro)=0, (25) 

where the parameter v has a meaning of the sound velocity in the center of the trap: 

^2^Ptf(%^A 
m m 

The substitution ttq = e^'^'^u{x) transforms (25) into the Legendre equation: 

.u{x) + ^( (l-^)^uix)] =0. (27) 



ffv"^ dx\\ R^J dx 

Since the Thomas-Fermi solution (23) is non-zero only at \x\ < R^ we shall consider (27) 
at X e [—Rc: Rc] C R, as well. After an analytical continuation u — > iE, equation (27) 
possesses the polynomial solutions, which are given by the Legendre polynomials Pn{x/Rc), 
if and only if 



7 



In other words, equation (27) leads to the spectrum of the low lying excitations: En = 

hVL\l ^ ) ^ 0, [18]. Notice that the corresponding equation for the homogeneous 

Bose gas is obtained after a formal limit 1/Rc in (27) at finite x. Provided the latter 
is still considered for the segment [—Rci -RJ 3 x with a periodic boundary condition for 
we arrive at the discrete spectrum of the following form: Ej^ — hvk, where k is the wave 
number, k — {Tr/Rc)n, n e Z. 



3 The two-point correlation functions 

Let us go over to our main task — to the calculation of the two-point thermal correlation 
function T{xi, Ti, X2, T2) of the spatially non-homogeneous Bose gas. We define it as the ratio 
of two functional integrals: 

Jij{xi,n)^{x,,r2)e^^^'^^VijV'ilj 
= J eSi^^^Wi^Vi. ' 

where the action S[ilj,ijj] is given by (2). 

We are interested in the behaviour of the correlators at the distances considerably smaller 
in comparison with the size of the domain occupied by the condensate. The main contribution 
to the behaviour of the correlation functions is due to the low lying excitations at sufficiently 
low temperatures [8, 11]. To calculate r(xi, Ti; ^2, T2), (29), we use the method of successive 
functional integration first over the high-energy excitations 'ijje, i^ei and then over the low- 
energy excitations ipo-, i/^o (see (4)). In the leading approximation, the correlator we are 
interested in looks, in terms of the density-phase variables, as follows [5,6]: 

r(xi,Ti;X2,T2) ~ J exp(^Ses[p,^]-i^{Xl,Ti)+iip{x2,T2) + 

+ i lnp(a;i,ri) + | In p{x2,T2)^VpV(fX (39) 
exp(5'efr[p, (p])VpV(p 



X 



-1 



where the integrand in the nominator is arranged in the form of a single exponential. Here 
Ses[p,^] is the effective action (21). 

Since the fiuctuations of the density are suppressed at sufficiently low temperatures [16], 
one can replace In p(xi, ri). In p{x2, T2) in (30) by In p^Fixi), In pTFix2), where pxF is defined 
by (23). In accordance with the variational principle suggested in [14], we estimate the 
functional integrals in (30) by the stationary phase method. Each of the integrals is charac- 
terized by its own stationarity point given by variation of the corresponding exponent. For 
the correlation function T{xi, ti, X2, T2), we obtain the following leading estimation: 

r(a;i,ri;a;2,T2) ~ \/ Pt:f{xi)pt:f{x2)^ 

X exp(-S'eff[po,¥'o]+'S'eff[pi,</2i] - 'npi{xi,Ti) + i</9i (a;2, T2)) , 

where the variables po; a-re defined by the extremum condition 5(5'eff[p, </?]) = 0, and 
therefore they just satisfy the Gross-Pitaevskii equations (22). The fields pi, </7i are defined 
by the extremum condition: 
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(s^s[p, (/?] - i(/7(xi, Ti) i</?(x2, T2)) = . (32) 
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The variational equation (32) leads to another couple of equations of the Gross-Pitaevskii 
type. One of these equations turns out to be a non-homogeneous equation with the S-hke 
source, while another one is a homogeneous equation. In fact, the homogeneous equation 
appears due to a requirement of vanishing of the coefficient at the variation 6p{x,r), while 
the non-homogeneous equation is defined by vanishing of the coefficient at the variation 

Sip{x,T). 

It can consistently be assumed that the solution pi{x, r) can be represented as a sum of 

PTpix) and of a weakly fluctuating part, provided the boundary is far from beginning of 
coordinates: pi{x,r) = Ptf{x) + 7ii{x,r). Therefore, the terms and dx^idx^i 

are small and can be omitted. Taking into account a linearization near the Thomas-Fermi 
solution, one can finally arrive at a couple of the following equations: 

^2 (33) 

-l^r'Kl^ d^{pT:Ydx^i) = i5{x - xi)5{t - Ti) - i5{x - X2)5{t - T2) . 

m 

Equations (33) lead [5,6] to the following equation for the variable </7i: 

-^^dl^pi + dx(pT:Y{.x)dxLpi) = i^|5(x-xi)5(T-ri) - 5{x - X2)5{t - T2)^, (34) 

where v means the sound velocity in the center of the trap (26), and prr is defined by (23). 
Now, with the help of (33) one can calculate the terms contributing into the exponent in 
(31), [5,6]: 

-5'eff[po, Vd] + 'S'eff[pi, <Pl\ ^ - ((/?i(xi, Ti) - (pi{x2, T2)) . (35) 

Substituting (35) into (31), one obtains the following approximate formula for the correlator: 

r(xi, Ti; X2, T2) ~ \/Ptf{Xi)ptf{x2) exp ^-^ (yJi(xi, Ti) - </7l(x2, T2))^ . (36) 

It is natural to represent solutions of equation (34) in terms of the solution G{x, r; x\ t') 
of the equation 



^ dlG{x, r; x\ r') + (^{l - ^) d^^, r; x', r')) = ^ 



6{x - x')6{t - t') . (37) 



Bearing in mind the homogeneous equation (25), we shall call (37) as non-homogeneous 
Stringari equation. As it is clear after [19], the Green function G(xi, ri; X2, T2) has a meaning 
of the correlation function of the phases: 

G{xi,Ti]X2, T2) = -((/?(xi, Ti)(/7(X2, T2)), (38) 

where the angle brackets in right-hand side should be understood as an averaging with 
respect to the weighted measure exp {Ses\p, </?]) VpVip. Using (38), it is possible to represent, 
eventually, the correlation function as follows [5,6]: 



r(a;i, Ti; X2, T2) ~ \/ p{xi)p{x2) exp {-\ {G{xi, n; X2, T2) + G{x2, T2; Xi, ri))) , (39) 

where /o(a;i), ]}{x2) are the renormalized densities. The solution G{xi^ti\X2^T2) of equation 
(37) is defined up to a purely imaginary additive constant, which has a meaning of a global 
phase. 
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The governing equations reported in [4] for the spatial dimensionahties d — 3,2, 1 are in a 
direct agreement with (33), provided the r-dependence is neglected in (33). The correlation 

functions of the phases are obtained in [4] without an influence of d^PTF dx^Pi as follows (the 
notation /(x, x') is used for them in [4], but G(x, x') is used below to keep contact with 
(38)): 

^''-'-''-- wp..(5)ix-x-r ^''=''>' 

^(^•'''' = 2^fidbs) i^-'^'i- ''"^^ 

where x, x' label spatial arguments at d = 3, 2, 1 and 5* = | (x + x'). In (40b), the thermal 
length At = h(3v is introduced, where v = \/~KJm is the sound velocity given by (26). It 
is already clear that the correlation functions can no longer depend on | x — x' | alone: 
they depend also on the center of mass coordinate S, consistent with the breakdown of 
translational invariance induced by the trap. 

Let us remind first the correlation functions in d = 3, 2. In this case, the points x = Xi 
and X = X2 in = i^i(x;xi,X2) are singular and introduce a divergence problem [8]. This 
difficulty can be avoided by considering a first-order coherence junction r^^^(xi,X2) which is 
defined in [4] as 

rW(x„x,)^— (41) 

where r(xi,X2) is given by (36). Then, r*^^^(xi,X2) is both finite and well defined because 
identically the same singularities appear [8] in a direct calculation of ('^o(xi)) and l^oi^'i))- 
We find that 

r(^)(xi,X2) - exp ( — ^zt4 7^ 1 ^ r J • (42) 

V47r/?^VpTF('5) I Xi - X2 17 



Evidently for d = 3, 



rW(xi,x2)^i+ ^ 



AnfJh^v'^PTFiS) I xi - X2 I 

for I xi — X2 |3> A/ {AttPH'^v^ptf^S)), thus indicating long-range order and long-range co- 
herence. The correlation length given by K/ Atk pii^v'^ pTp{S)) is therefore a slowly varying 
function of Xi and X2. Notice that we have assumed Xi and X2 are not close to the bound- 
aries. 

In the case d = 2, the correlations decay by a power law for arbitrary small temperatures, 
rW(x„X2)c.(^-^) . (43) 

\ I Xi — X2 I / 

The exponent of this power-law is proportional to T so that, at very low temperatures, 
correlations may thus prevail over almost macroscopic distances. At T = exactly, T^^^ will 
have a non-zero off-set due to the presence of a true condensate. 
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For d — 1, one can obtain for G instead of (40c) an exact expression as follows (i.e., the 
term O^Ptf d^^Pi in the governing equations is now accounted for): 



^(x, x') = G{x, x') 



gRc 



p{2hv) 



In 



{l + \x- x'\/{2R,)Y -{x + .T')V(4i?2) 
{l-\x- x'\/{2Rc)Y - (x + x')V(4i?2) 



(44) 



Then, we obtain r(xi,X2) = r(xi, X2) in the following form: 



1 + |a;i - X2\/Rc - X1X2/RI 



1 - \xi - X2\/Rc — X1X2/RI 

In the limit | Xi — X2 |<^ {Rc, S}, we obtain r(xi, X2) in the form: 

A 



(45) 



V{xi,X2) ~ \/ Ptf{xi)ptf{x2) exp 



2j3h^v^p^Y{S) 



\Xi - X2 



(46) 



The correlation length depends now on Xi and X2. Without the trap, the correlation length 
reduces to 2/3/i^pTF(0)/m, where Ptf(O) is the density of the ground state of a homogeneous 
system at zero temperature, in complete agreement with the exact solution [3]. 

It is obvious that the correlation functions (43) and (46) vanish for large separation of 
the arguments at arbitrary small temperatures T > 0, and that there is no long-range order 
in (i = 2 or (i = 1. The correlation functions (42), (43), (46) coincide with those obtained 
under translational invariance without the trap to the extent that there is now an additional 
factor ptf(>S') in the exponents. True long-range order arises only in d = 3. 



4 The asymptotics of the correlation functions 

Therefore, the problem concerning the study of the asymptotical behaviour of the two-point 
thermal correlation function r{xi,ri;x2,T2) given by the representation (39), is reduced to 
solution of the non- homogeneous Stringari equation (37). The corresponding answer (or its 
asymptotics) should be subsequently substituted into (39). In the present section, we shall 
obtain explicitly solutions of (37), and we shall consider the corresponding asymptotics of 
r{xi, Ti, X2, T2). Let us begin with the limiting case of a homogeneous Bose gas. 

4.1 The homogeneous Bose gas 

The homogeneous case is given by V{x) = 0, and the related equation appears from (37) at 
1/Rc ^ 0: 

^ d^,G{x, r; x', r') + dlG{x, r; x', r') ^-^5{x- x')5{r - r') . (47) 



We consider (47) for the domain [—Rc,Rc] x [0,(3] 3 {x,t) with the periodic boundary 
conditions for each variable. The 5-functions in right-hand side of (47) are treated as the 
periodic 5-functions. This allows us to represent the solution of this equation as the formal 
double Fourier series: 

g'-'-^-''-'>-(2^)^ + ' 

oj.k 
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where uj — {2n/P)l, I e Z. The notation for the energy Ek = hvk, where k = {n/Rc)n, 
n e Z, is used in (48). Besides, the representation (48) requires a regularization, which 
consists in neglect of the term given by a; = A; = 0. 

Let us deduce from (48) two important asymptotical representations for the Green func- 
tion. Then, in the limit of zero temperature and of infinite size of the domain occupied by 
the Bose gas, one can go over to the asymptotics of r(a;i, ri; X2, T2). When a strong inequality 
= ksT » hv/Rc is vahd, we obtain: 



G{x, r; x', t') 



9 



2nhv 



\n<2 



sinh — — (\x — x'\ + ihv{T — t')) 
npv ^ ' 



9 



\x — X 



l\2 



(49) 



where — a;'| < 2i?c, |r — r'| < /?, and C is some constant, which is not written explicitly. 
When an opposite inequality = ksT <^ hv/Rc is valid, we obtain: 



G{x, r; x' , t') ~ 



ln<^ 2 



2TThv 



ITT 

sinh — — (\x — x'\ + ihv{T — r')) 



r'|2 



4/9i?c 



+ C', 



(50) 



where \x — x'\ < 2R,., |r — r'| < /5, and C is another constant. 

Let us substitute the estimate (49) into the representation (39) and take simultaneously 
the limit (5hv/ Rc (the size is growing faster than inverse temperature). Then, we obtain 
the following expression for the correlator in question: 



r(xi,Ti;x2,T2) ~ ^/p{xi)p{x2) sinh^^dxi -a;2| +i^v(Ti -T2)) 



-g/2'Khv 



(51) 



Further, applying the relation (50) and taking the limit R^./ {j3hv) — >■ (the inverse temper- 
ature grows faster than the size), we obtain for r(a;i, ri; 2:2, T2): 



T{xi,Ti,X2,T2) ~ sj p{xi)p{x2) 



sinh 



ITT 
2Rr 



\Xi - X2I + \hv{Ti - T2)) 



-gj2-Khv 



(52) 



It follows from (51) and (52), that in the limit of zero temperature, {h^v) — > 0, and of 
infinite size, 1/Rc ^ 0, the two-point correlation function behaves like 



r(xi,Ti;X2,T2) 



^/p{xl)p{x2] 



\xi - X2\+ ihv(ri - T2 



(53) 



The latter formula is valid in the limit Phv/Rc 0, as well as in the limit Rc/{(3Hv) — > 0. 
In (53) 6 denotes the critical exponent: 6 = 2iThv/g, and the arguments Xi and X2, Ti and 
T2 are assumed to be sufficiently close each to other. Using the notations v = kjm for the 
sound velocity and p = kj g ioY the density of the homogeneous Bose gas, we obtain for the 
critical exponent the following universal expression: 



2'nhp 



mv 



(54) 
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4.2 The trapped Bose gas. High temperature case: ksT ^ hv/Rc 

Let us turn to non-homogeneous Bose gas described by equations (l)-(3). Now, we should 
consider the non-homogeneous Stringari equation (37) for the arguments {x, r) e [—Rc, Rc] x 
[0, /3] with the periodic boundary condition only respectively to r (contrary to equation (47), 
S{x — x') is a usual Dirac's 5-function supported at the point e R). The Green function 
satisfying (37) can be written as a formal Fourier series: 

G{x, r; x\ r') = - Y, e^^^^^-^'^G^^x, x') , (55) 

where uj = (27r//3)/, / e Z. The spectral density Gt^{x,x') in (55) is then governed by the 
equation 



+ ^ ( (l - g A OM, .')) = ^ . (56) 



Solution of equation (56) can be obtained in terms of the Legendre functions of the first 
and second kind, Py{x/Rc) and Qi,{x/Rc), [20], which are linearly independent solutions of 
the homogeneous Legendre equation (27). As a result we get [5,6]: 



G^{x, x') = ^G^{x, x') + i^G^ix, x') , (57) 



where 



^ ||j.(.-.'){a(f)P.(|;)-a(|^)P.(| 



(68) 



V looks as follows: 



2 V4 Khv, 



and e{x — x') is the sign function e{x) = sign(a;). Validity of the solution (57), (58) can 
be verified by direct substitution into (56), where an expression for the Wronskian of two 
linearly independent solutions P^, and Q^, [21] should be used. 

The Green function (57) can be represented in the form, which allows to study the 
corresponding asymptotical behaviour. When the coordinates Xi, X2 arc chosen to be far 
from the boundary of the trap, Xi, X2 -C Rc, but at the same time the inequalities |a;i — a;2| <^ 
I {xi + X2) and \xi — X2\ <^ Rc are valid, the corresponding limit should be called as quasi- 
homogeneous. In the case of strong inequality j3~^ = ksT ^ hv/Rc, we approximately obtain 
for non-zero frequencies: \uj\ ^ hv/{2Rc). Using the standard asymptotics of the Legendre 
functions [20,22], we determine the behaviour of Gi^{x,x') in the quasi-homogeneous limit 
at large |a;| as follows: 

n f '\ ^ exp{-{hv)-^\uj\\x - x'\) 

G^{x,x):^-— — ^ — 59 

Here S means a half-sum of the spatial arguments of the correlator, S = ^ {xi + X2), and v 
is given by (26). 
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Further, we find that the term P~^Gq{x, x') in (55) is just given (in the quasi-homogeneous 
hmit) by G(x, x') (40c). Using (40c) and (59) for evaluation of the series (55), one obtains 
the answer: 



Therefore, the Green function (39) takes the following form at j3~^ 3> hv/Rc- 



(60) 



^ Vp(xi)p(x2) , . 

r(a;i,ri;a;2,r2) ~ ^ , (61) 

sinh — - (\xi - X2\ + ihv{Ti - T2)) 
where the critical exponent ^(5") depends now only on the half-sum of the coordinates S: 

mv 

The result (61), which is valid for the spatially non- homogeneous case, is in a correspon- 
dence with the estimation (51) obtained above for the homogeneous Bose gas. Therefore, 
the expression (61) is also concerned with vahdity of the condition that the size of the do- 
main occTipicd by the Bose condensate grows faster than inverse temperature, i.e., with the 
condition hj3v/ Rc 0. 

The relation (61) can be simplified for two important limiting cases. Provided the con- 
dition 

npv npv 

is fulfilled in the quasi-homogeneous case, we obtain from (61) that the correlator decays 
exponentially: 



r(a;i,ri;a;2,r2) ~ \/p{xi)p{x2) exp ~ ^^1 + ihv{Ti - T2)\^ , 



(64) 



2(3h?v'^PTF{S) 



The correlation length ^,{S) is defined by the relation (64), which depends now on the half- 
sum of the coordinates: 

5(s).5^(,(s) = ?5^*5i(£). (65) 

TT m 

In an opposite limit, 

\Xl-X2\ |n-T2| Rc 

' R ^ ^ ^ ^ ' (66) 

npv p npv 

the asymptotics of T{xi,ti;x2,T2) takes the following form: 

r(a;i, n; X2, T2) ~ ^— . (67) 

|F1 - ^2] +mv{Ti - T2)\ 

The obtained asymptotics (67) is analogous to the estimation (53), which characterizes the 
spatially homogeneous Bose gas. But the critical exponent 6{S) in (67) differs from 6 (54), 
since the latter does not depend on the spatial coordinates. 
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4.3 The trapped Bose gas. Low temperature case: fc^T <C fw / Rc 



Let us go over to another case, which also admits investigation of the asymptotical behaviour 
of the two-point correlator r(a;i, ri; 0:2, T2). We begin with the non- homogeneous Stringari 
equation (37), which can be rewritten in the following form: 



dlG{x, r; x' , t') + (l - ^) d{.x/Rc)G{x, r; x' , t') 



9_ 

Rr 



(68) 



where the notation a = Rc/ {hv) is introduced. The asymptotical behaviour can be inves- 
tigated in two cases: P/a <^ 1 (the previous subsection) and P/a 1 (see below). The 
functions 

X 
-Rr 



n>0, 



where Pn{x/Rc) are the Legendre polynomials, constitute a complete orthonormal system in 
the space L2 [—Rc, Rc]- This fact allows to obtain the following representation for the Green 
function G{x, r; x', r') in the form of the generalized double Fourier series: 



G(x,x;x'./)^(^ 



00 ,1 



oj n=Q 



UJ- 



x 

Rr 



P 



(69) 



Summation ^ in (69), as well as in (48) and (55), goes over the Bose frequencies, and the 
following notation for the energy levels (28) is adopted: 



n{n + l) _ v/ri(riTT) 



a 



(70) 



After summation over the frequencies and after regularization consisting in neglect of the 
term corresponding to zero values of cu and n, G{x, r; x', t') (69) takes the form: 



P X 2 ~ COs(^/) 



+ 



00 I 1 
n+ i 



n=l 



Ern 



Pn{^)Pn{^) (coth(l (5E^) cosh(E„AT) - sinh(£;„Ar)) 



(71) 



where Ar = |r — r'|. The representation (71) can be studied in both cases: l3/a <^ 1 and 
13 /a 3> 1. For instance, using (71) at coinciding arguments r = r' to obtain T{xi^t\X2,t) 
in the limit 13 /a -C 1, we just obtain [6] the representations (45) or (under the quasi- 
homogeneity condition) (46). 

Let us turn to the case ^ 1, where PEn S> 1, Vn > 1. In other words, let us suppose 
that UbT <^ En and, so, ksT <^ Kl. Then, one obtains from (71): 



G{x, r; x', r') 



4i?. 



(I- 



Ar\2 1 
~p) ~12 



(72) 
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Notice that a difference between two neighbouring energy levels (70) can be estimated. 
After some appropriate series expansions, which are valid at Vn > 1, one obtains: 



En+l — E„ 



1 + 



+ ... 



(73) 



Equation (73) demonstrates that the levels (70) are approximately equidistant provided the 
inverse powers of n are neglected in (73) for sufficiently large n > uq. In its turn, the 
following estimation is also valid: 



n + 1/2 
\A(^H^ 



1 + 



+ 



(74) 



It is remarkable that the terms oc are absent both in (73) and (74). Let us remind that 
leading asymptotical estimations, which are obtainable with so-called logarithmic accuracy, 
are important for physical applications. The problem at hands just admits an estimation 
with leading logarithmic accuracy, since inverse powers of n can be omitted with the same 
(and good) accuracy in (73) and (74) at sufficiently large n. Convergency of the series (72) 
is not affected in this situation, since Ar is non-zero. 

It is known that at sufficiently large n, the following asymptotics for the Legendre poly- 
nomials Pn is valid [21]: 



Pn(cOS^?) 



TTO sin 1? 



cos 



(n + |)i?- +0(n-=^/'), 0<i?<7r. 



(75) 



n=no 



Let us split the sum over n in (72) into two parts: and . Further, let us assume 

n=l n=no+l 

that riQ is large enough to substitute at n > no the Legendre polynomials Pn by their 
asymptotical expressions given by (75) (and denoted below as P„ (cos ■»?)). Then, using (73) 
and (74), we can put G{x,t]x\t') (72) into the following approximate form [6]: 



G{x, r; x', t') 

X exp {^—■\/n{n -\- 1) 



P J 12 



9 
2hv 



no 

E 

n=l 



a ) 



\/n(rrTT) 

-^»(^)^»(i)'='K-(''+^) 



X 
Rr 



a J 



(76) 



2hv 



oo , 



n=l 



where t = exp(— Ar/a), and Hq is the number, which is fixed (its specific value is forbidden 
to go to infinity). Expression (76) is valid when r and r' arc close either to zero or to (3. 
Besides, we assume that r 7^ r' in order to keep convergency of (76), while the term omitted 
can be estimated [6]. 

Provided a smallness of x/Rc, x' j R^ and At /a is taken into account, a leading logarithmic 
behaviour of the scries in (76) can be established by the standard tools [21]. Since in the 
logarithmic approximation the first two terms in (76) are less important in comparison to 
the third one, we write down the leading contribution to the Green function G(x, t;x',t') 
in the quasi-homogeneous limit as follows: 



G{x, r; x\ t') 



-A 



In 



Rc 



2'Khv Ptf{S) \x — x'\ + ihv{T — t') 



(77) 
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where pxF is given by (23), 5" is a half-sum of x and x', and it is assumed that 



I \x — x'\ + ihviT — t') I 
= ^ ^—-^ ^ < 1 . (78) 

Besides, at sufficiently large no, in our consideration it is more appropriate to keep n instead 
of n + I in (75). Expression (77) takes place provided the following conditions of validity of 
the logarithmic estimation are respected: 

111 

Kno< — < — In — . (79) 

XL-jfi 1/'s(( 

A specific value of Uq can be related to the size of the trap 7?^: at a restricted range of 
deviations between the spatial arguments x and x\ increasing of Rc implies increasing of 
an upper bound for admissible values of no. However, due to (79), the estimation obtained 
for G{x,t]x\t') (77) does not depend explicitly on a specific choice of no. In the limit 
1/Rc — > 0, the total coefficient in front of the logarithm in (77) acquires the value —1/9, 
where the critical exponent 9 is defined like in (53), (54). 

Eventually, we obtain the following estimation for the two-point correlator 
r(xi,ri;x2,r2): 

j.( X v^p(a;i)p(a;2) , . 

r(xi, n; X2, T2) ~ ■ -rneis) ' (^o) 

|Fi - X2I ihv{Ti - T2)\ 

where the notation for the critical exponent ^(5*) is given by (62). 

The estimation obtained (80), where the critical exponent is 9{S) (62), constitutes the 
main result of the present subsection devoted to the case given by ksT <^ hv/Rc- Prom a 
comparison with the spatially homogeneous Bose gas, one can see that now the derivation 
of the estimate (80) is just analogous to a transition from the relation (52) to the final 
asymptotics (53). Then, validity of the corresponding limiting condition Rc/{hj3v) 
means that the result (80) is also due to the fact that the condensate boundary Rc increases 
slower than the inverse temperature. 

Therefore, under the different conditions, ksT ^ hv/Rc and ksT <C hv/Rc, we demon- 
strated in the present section that the behaviour of the correlator in the limit of zero tem- 
perature, {h(3v)~^ 0, and of infinite size of the trap, 0, is given by the coinciding 
estimations (67) and (80), i.e. the two-point correlation function T{xi,Ti;x2, T2) has a unique 
power-law behaviour in this limit. 



5 Conclusion 

The model considered describes a spatially non-homogeneous weakly repulsive Bose gas 
subjected to an external harmonic potential. The functional integral representation for the 
two-point correlation function is estimated by means of the stationary phase approximation. 
The main results are obtained for the case when the size of the domain occupied by the 
quasi-condensate increases, while the temperature of the system goes to zero. In the one- 
dimensional case, the behaviour of the two-point correlation function at zero temperature is 
governed by a power law. However, in contrast with the case of spatial homogeneity of the 
Bose gas, the corresponding critical exponent depends on the same spatial arguments as the 
correlator itself. It is just the presence of the external potential which is responsible for the 
non-homogeneity of the critical exponent. 
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